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Abstract 

In global seismology Earth's properties of fractal nature occur. Zyg- 
mund classes appear as the most appropriate and systematic way to mea- 
sure this local fractality. For the purpose of seismic wave propagation, we 
model the Earth's properties as Colombeau generalized functions. In one 
spatial dimension, we have a precise characterization of Zygmund regular- 
ity in Colombeau algebras. This is made possible via a relation between 
moUifiers and wavelets. 



1 Introduction 

Wave propagation in highly irregular media. In global seismology, (hyperbolic) 
partial differential equations the coefficients of which have to be considered 
generalized functions; in addition, the source mechanisms in such application 
are highly singular in nature. The coefficients model the (elastic) properties of 
the Earth, and their singularity structure arises from geological and physical 
processes. These processes are believed to reflect themselves in a multi-fractal 
behavior of the Earth's properties. Zygmund classes appear as the most appro- 
priate and systematic way to measure this local fractality (cf. ||, Chap.4]). 

The modelling process and Colombeau algebras. In the seismic transmission 
problem, the diagonalization of the first order system of partial differential 
equations and the transformation to the second order wave equation requires 
differentiation of the coefficients. Therefore, highly discontinuous coefficients 
will appear naturally although the original model medium varies continuously. 
However, embedding the fractal coefficient first into the Colombeau algebra en- 
sures the equivalence after transformation and yields unique solvability if the 
regularization scaling 7 is chosen appropriately (cf . [0, ^ ) . We use the frame- 



work and notation (in particular, Q for the algebra and An for the moUifier sets) 
of Colombeau algebras as presented in ||T^ . 

An interesting aspect of the use of Colombeau theory in wave propagation is 
that it leads to a natural control over and understanding of 'scale'. In this 
paper, we focus on this modelling process. 

2 Basic definitions and constructions 

2.1 Review of Zygmund spaces 

We briefly review homogeneous and inhomogcncous Zygmund spaces, CKR™) 
and C;^(R™), via a characterization in pseudodifferential operator style which 
follows essentially the presentation in Sect. 8.6. Alternatively, for practi- 
cal and implementation issues one may prefer the characterization via growth 
properties of the discrete wavelet transform using orthonormal wavelets (cf. [^). 

Classically, the Zygmund spaces were defined as extension of Holder spaces 
by boundedness properties of difference quotients. Within the systematic and 
unified approach of Tricbcl (cf. jTs], ^) we can simply identify the Zygmund 
spaces in a scale of inhomogeneous and homogeneous (quasi) Banach spaces, B^^ 
and B;^ (s e M, < p, g < oo), by C|(M™) = S^«,(M™) and (7|(R™) = B^^. 
Both C^^(M") and CKM™) are Banach spaces. 

To emphasize the close relation with moUifiers we describe a characterization of 
Zygmund spaces in pseudodifferential operator style in more detail. 

Let < a < 6 and choose tpo £ T>{M.), ipo symmetric and positive, (po{t) = 1 
if \t\ < a, (po{t) = if \t\ > h, and Lp^ strictly decreasing in the interval {a,b). 
Putting (p{^) = tpodCI) for ^ e M™ then defines a function tp e 2?(R"). Finally 
we set 

V'(0 = -(e|grad^(0) 

and note that if a < |C| < 6 then i/>(^) = -iPo{\^\)\£.\ > 0. We denote by X(M'") 
the set of all pairs {(p,4') G I?(R'")^ that are constructed as above (we usually 
suppress the dependence of on a and b in the notation). 

We are now in aposition to state the characterization theorem for the inhomo- 
geneous Zygmund spaces as subspaces of 5'(M'"). It follows from |Q, Sec. 2.3, 
Thm. 3 or, alternatively, from ||], Sec. 8.6. Note that all appearing pseudod- 
ifferential operators in the following have x-independent symbols and are thus 
given simply by convolutions. 

Theorem 1. Assume that a < 1/4 and 6 > 4 and choose {(p,tp) G M{M."^) 
arbitrary. Let s G R then u G iS'(R™) belongs to the inhomogeneous Zygmund 
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space of order s CJ(R™) if and only if 

(1) \u\cs MD)u\\^^ + sup (t-'UitD)u\\^^) < oo. 

o<t<i ^ ^ 

(Note that we made use of the modification for (7 = cxo in [|^, equ. (82).) 

Remark 2. (i) \u\c= defines an equivalent norm on C^. In fact that all 
norms defined as above by some {tpjtjj) £ A^(R™) are equivalent can be 
seen as in Lemma 8.6.5. 

(ii) If s e E+ \ N then C|(IR'") is the classical Holder space of regularity 
s. Denoting by [sj the greatest integer less than s it consists of all [sj 
times continuously differentiable functions / such that d^f is bounded 
when |a| < [sJ and globally Holder continuous with exponent s — [sJ if 
\a\^[s\. 

(iii) Due to the term ||(^(£')u||loo the norm \u\c^ is not homogeneous with 
respect to a scale change in the argument of u. 

(iv) If u e L°°(M") then (cf. @, Sect. 8.6) 

00 

ij){D /t)u{x)— for almost all x. 

1 

Using ip{(,) = tp{^/t)/tdt this can be rewritten in the form u{x) = 
^p{D/t)u{x)/t dt and resembles Calderon's classical identity in terms 
of a continuous wavelet transform (cf. [D, Ch. 1, (5.9) and (5.10)). 

(v) In a similar way one can characterize the homogeneous Zygmund spaces as 
subspaces of 5'(M'") modulo the polynomials V. A proof can be found in 
Q, Sec. 3.1, Thm. 1. We may identify S' /V with the dual space 5^(R'") 
of 5o(IR") = {/ e 5(R™) I (9"/(0) = OVa G N^'}, the Schwartz functions 
with vanishing moments, by mapping the class u + "P with representative 
u <E S' to u \sg. Assume that a < 1/4 and b > 4 and choose ip G ViW^) 
as constructed above and let s G M and u G iS'(M™). Then u\so belongs 
to the homogeneous Zygmund space C'*(IR™) of order s if and only if 

(3) |u|c= := sup (t-''\\tP{tD)u\\^^)<oo. 

0<t<oo ^ ■' 

(Note that we use the modification for g = 00 in |l2j, equ. (16).) 
2.2 The continuous wavelet transform 

Following H we call a function g G L^iW) D L°°(M'") with / g = a wavelet. 
We shall say that it is a wavelet of order k {k E Nq) if the moments up to order 
k vanish, i.e., J x"g(x)dx = for \a\ < k. 
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The (continuous) wavelet transform is defined for / G Lp(M™) (1 < p < oo) by 
(e>0) 

(4) WJ{x,s)^ I f{y)l--g{y—^)dy^f*{lUx) 

where we have used the notation g{y) = g{—y) and gdy) = g{y/e)/e'^. By 
Young's inequahty Wgf{.,e) is in 17 for all e > and Wg defines a continuous 
operator on this space for each e. 

li g & Cc(K™) we can define Wgf for / G Lj'Q^(R'") directly by the same formula 
(|). If 5 e 5o(M") then Wg can be extended to as the adjoint of the 

wavelet synthesis (cf. [0, Ch. 1, Sects. 24, 25, and 30) or directly by S'-S- 
convolution in formula 

Remark 3. If / is a polynomial and g E So it is easy to see that Wgf — 0. 

In fact, /, g, and Wgf are in S' and (Wg/(.,e)) ~ fg{e.). Since g is in So 
the Fourier transform g{e.) is smooth and vanishes of infinite order at 0. But 
/ has to be a linear combination of derivatives of 5o implying fg{e.) = 0. 
Therefore the wavelet transform 'is blind to polynomial parts' of the analyzed 
function (or distribution) /. In terms of geophysical modelling this means that 
a polynomially varying background medium is filtered out automatically. 



2.3 Wavelets from mollifiers 

The Zygmund class characterization in Theorem ^ (and remark ^(v)) used 
asymptotic estimates of scaled smoothings of the distribution which resembles 
typical moUifier constructions in Colombeau theory. In this subsection we relate 
this in turn directly to the wavelet transform obtaining the well-known wavelet 
characterization of Zygmund spaces. 

Let X G 5(18.™) with / x = 1 and define the function by 



(5) jl{x) := mx{x) + (.T|grad x(a;)) . 

Then /i is in S{W^) and is a wavelet since a simple integration by parts shows 
that 



(-1)1"! / ^l{x)x'^dx = j fi{x){-x)°' dx 

P m „ 

= ,-l)H«/.",,.,^i:o, = (-l)-'|o|/."xW.. 



i=i 



J fi — and if \a\ > we have J x" ^{x) dx = if and only if J x"x(x) dx ~ 0. 
Therefore /i defined by (|^) is a wavelet of order N if and only if the moUifier x 
has vanishing moments of order 1 < la] < N. 



4 



Furthermore, by straightforward computation, we have 

(6) (A)e(x) = -£^(Xe(a:)) 

yielding an alternative of in the form fl{x) — ~-j^{Xe{^)) \s=i- 

If {(p,ip) G A1(R™) arbitrary and x, M are the unique Schwartz functions such 
that X = 'P and — ^p, then straightforward computation shows that /z and x 
satisfy the relation (||). Therefore since n is then a real valued and even wavelet 
we have for u £ S' 

i:(tD)u{x) = f^U-) * u{x) = W^,u{x, t) . 

Hence the distributions u in the Zygmund class C|(M™) can be characterized 
in terms of a wavelet transform and a smoothing pseudodifferential operator by 
< oo and supo<t<i {t~^\\W^u{.,t)\\^^ < oo. We have shown 

Theorem 4. Let (x, m) G A^(M")- A distribution u e iS'(M") belongs to the 
Zygmund class C| (M™ ) if and only if 

(7) ||u*xIIl-<oo and \\W^u{.,r)\\^^ ^ 0{r') {r ^ 0). 

Remark 5. (i) Observe that the condition on x implies that x and hence 
fi can never have compact support. If this characterization is to be used 
in a theory of Zygmund regularity detection within Colombeau algebras 
one has to allow for moUifiers of this kind in the corresponding embedding 
procedures. This is the issue of the following subsection. Nevertheless 
we note here that according to remarks in ^, (2.2) and (3.1)] and, more 
precisely, in ||^, Ch.3] the restrictions on the wavelet itself in a charac- 
terization of type (|^) may be considerably relaxed — depending on the 
generality one wishes to allow for the analyzed distribution u. However, 
in case m = I and u a function a flexible and direct characterization (due 
to Holschneider and Tchamitchian) can be found in 0], Sect. 2.9, or |H, 
Sect. 4.2. 

(ii) There are more refined results in the spirit of the above theorem describing 
local Holder (Zygmund) regularity by growth properties of the wavelet 
transform (cf. in particular ||^, Sect. 4.2, and |^). 

(iii) The counterpart of (|^) for Lj^^-functions in terms of (discrete) multireso- 
lution approximations is Sect. 6.4, Thm. 5. 

3 Colombeau modelling and wavelet transform 
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3.1 Embedding of temperate distributions 



We consider a variant of the Colombeau embedding : I?'(R'") Q{W'^) 
that was discussed in [Q, subsect. 3.2. As indicated in remark ||,(i) we need to 
allow for moUifiers with noncompact support in order to gain the flexibility of 
using wavelet-type arguments for the extraction of regularity properties from 
asymptotic estimates. On the side of the embedded distributions this forces 
us to restrict to iS', a space still large enough for the geophysically motivated 
coefhcients in model PDEs. 

Recall (Q, Def. 11) that an admissible scaling is defined to be a continuous 
function 7 : (0,1) R+ such that 7(r) — 0(l/r), 7(r) 00, and 7(sr) = 
0(7(r)) if < s < 1 (fixed) as r ^ 0. 

Definition 6. Let 7 be an admissible scaling, x G 5(M™) with Jx = I7 then 
we define q : 5'(]R") ^(K'") by 

(8) Ll{u) = cl[(u * x^'i^, " e 5'(M'") 

where 

{9) x'^{<l>,x) ^-f{l{(f>Q)rx{l{KM)x) if (I) = (f>o (1)0 with 00 e A- 

q is well-defined since (i^, x) — > u*x'^((/), x) is clearly moderate and negligibility 
is preserved under this scaled convolution. By abuse of notation we will write 
q(u){4', x) for the standard representative of tx(")- 

The following statements describe properties of such a modelling procedure re- 
sembling the original properties used by M. Oberguggenberger in [0, Prop. 1.5, 
to ensure unique solvability of symmetric hyperbolic systems of PDEs (cf. |^ , 
. The definition of Colombeau functions of logarithmic and bounded type is 
given in ||ll|, Def. 19.2, the variation used below is an obvious extension. 

Proposition 7. 

(i) q : 5'(M'") Q(U."'') is linear, injective, and commutes with partial 
derivatives. 

(ii) Vm e 5'(K"): q{u) « u. 

(iii) If M e Ty^^(K™) then q{u) is of ^-type, i.e., there is iV e No such that 
for aU (j) e ^Ar(R") there exist C > and 1 > 77 > 0: 

(10) sup \iliu){^,,y)\<Nj{Ce) < e < 77. 

(iv) If u G L°°(]R™) then ^{u) is of bounded type and its first order derivatives 
are of 7-type. 



6 



Proof, ad Is clear from Xe •= X^i't'ei-) 5 in S' as e ^ Q and the 

convolution formula. 

ad (in): Although this involves only marginal changes in the proof of pO| , 
Prop. 1.5(i), we recall it here to make the presentation more self-contained. 

Let u — uq + Sjli '^i^j with Uj e L°° {j — 0, . . . ,m) then with 7^ := 7(eZ(0o)) 

m 

\U * Xe{x)\ < \\uo * XeIIl- + Yl * 

m 

< ho|lL~llxllLi+7s5Zll^jllL-ll^jXllLi 

where the expression within brackets on the r.h.s. is bounded by some constant 
M, dependent on u and x only but independent of (f), as soon as e < 77 with 
r] chosen appropriately (and dependent on AI, u, x, and 0). Therefore the 
assertion is proved by putting N > M and C = l{4>o). 

ad (iv): is proved by similar reasoning □ 

In particular, we can model a fairly large class of distributions as Colombeau 
functions of logarithmic growth (or log- type) thereby ensuring unique solvability 
of hyperbolic PDEs incorporating such as coefficients. 

Corollary 8. 

(i) If 7(e) = log(l/e) then tJ(VF-i'°°) C {[/ e | [/ is of log-type } and 
lI{L°°) C {U e g \ U of bounded type and d"U of log-type for |a| = 1}. 

(ii) If u e W^'''°°(R"^) for A; G Nq then ij(u) is of 7''-type. In particular, there 
is an admissible scaling 7 such that l^{u) and all first order derivatives 
dji^{u) {j — 1,. . . , to) are of log-type. 

3.2 Wave front sets under the embedding 

One of the most important properties of the embedding procedure introduced in 
was its faithfulness with respect to the microlocal properties if 'appropriately 
measured' in terms of the set of 7-regular Colombeau functions Q'^iW"'^) (Q, 
Def. 11). But there the proof of this microlocal invariance property heavily used 
the compact support property of the standard moUifier x which is no longer true 
in the current situation. In this subsection we show how to extend the invariance 
result to the new embedding procedure defined above. 
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Theorem 9. Let w G 5'(R'"), 7 an admissible scaling, and x G <S(]R'") with 
/ X = 1 then 

(11) WFJiiliw)) = WF{w). 

Proof. The necessary changes in the proof of Thm. 15, are minimal once we 
established the following 

Lemma 10. If e P(E'") and v S with supp((p) n supp(i;) = then 

Proof. Using the short-hand notation Xe = X'^i4'ei ■) and 7^ = 7(£Z((/)o)) we have 

Hence we need to estimate terms of the form 7!"' {v, d^xileix — •))) when x G 
supp((y5) —: K. Let 5 be a closed set satisfying supp(ti) C 5* C M™ \ K and 
put d = dist(S', K) > 0. Since w is a temperate distribution there is G N and 
C > such that 

7l"l|(t;,9"x(7.(^--)))l <C7l"l E sup|9'^(a"x(7e(a:-y)))|. 

X G iS implies that each term in the sum on the right-hand side can be estimated 
for arbitrary fc G N by 

sup |5'^+"x(7.(^ - y))) l7i"' < 4"^ supCfe(l + 7.N - yj)-^ < cWf^-'^/d'' 
yes yes 

if X varies in K. Since \a\ + \a\ < \/3\ + N we obtain 

with a constant C depending on fc, v, ip, d, and x but k still arbitrary. Choosing 
fc = for example, we conclude that </? ■ iJJ(i') has a uniform 7j-growth over 
all orders of derivatives. Hence it is a 7e-regular Colombeau function. □ 

Referring to the proof (and the notation) of Q, Thm. 15, we may now finish 
the proof of the theorem simply by carrying out the following slight changes in 
the two steps of that proof. 

Ad step 1: Choose ip V such that = 1 in a neighborhood of supp(iy9) and 
write 

ip{w * Xe) = 'fiiitpw) * Xe) + "^(((1 - 1p)w) * Xe ■ 

The first term on the right can be estimated by the same methods as in § and 
the second term is 7-regular by the lemma above. 
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Ad step 2: Rewrite 

ipw = ipi/jW = ip{%ljW — {ipw) * Xe) + * Xe) 

and observe that the reasoning of Q is apphcable since Y^'^{ip i'^{'4jw)) C Y.'^{ip l^{w)) 
by the above lemma. 

□ 

3.3 The modelling procedure and wavelet transforms 

Simple wavelet-mollifier correspondences as in subsection 2.3 allow us to rewrite 
the Colombeau modelling procedure and hence prepare for the detection of orig- 
inal Zygmund regularity in terms of growth properties in the scaling parameters. 

A first version describes directly but involves an additional nonhomogeneous 
term. 

Lemma 11. If x G iS(K'") has the properties Jx = 1 a-nd / x°'x{^)dx = 
(0 < I a I < N) then /I — ~-§^{Xe) \e=i defines a wavelet of order N and we have 
for any / e 5'(]R") 

1 

/■ dr 

(12) il{m,x)^ f*x{x)+ J WJix,r)-. 

1/7(K0)) 

Proof. Let e > then eq. (|) implies Wf,f{x,e) = f*{fi)e{x) = -e£{f*Xe{x)) 
and integration with respect to e from 1/7(^(0)) to 1 yields 

1 

/■ de 
- J WJ{x,e)-^f*x{x)~q{f){<l>,x). 

l/7(iW)) 

□ 

A more direct mollifier wavelet correspondence is possible via derivatives of 
instead. 

Lemma 12. If x G iS(M™) with / X = 1 then for any a e Nq with |a| > 

(13) Xo.ix) = (d^xYix) 

is a wavelet of order |a| — 1 and for any / G iS'(M™) we have 

(14) d"q{f){^,x) = 7(^(0))l"l W^Jix, — 1-) . 
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Proof. Let |/3| < |a| then / xl^D°'x{x)dx = (--D)'^(rx(C)) l«=o= which 
proves the first assertion. The second assertion follows from 



d'-iUm) = 7'"'+™/ * a"x(7.) = 7'"'/ * (7™(5"x)"(7-))~ 
with the short-hand notation 7 = 7(2(0)). □ 

Both lemmas ^ and |l2| may be used to translate (global) Zygmund regularity 
of the modeled (embedded) distribution / via Thm. |7| into asymptotic growth 
properties with respect to the regularization parameter. To what extent this 
can be utilized to develop a faithful and completely intrinsic Zygmund regularity 
theory of Colombeau functions may be subject of future research. 



4 Zygmund regularity of Colombeau functions: 
the one-dimensional case 

If we combine the basic ideas of the Zygmund class characterization in 2.3 with 
the simple observations in 3.3 we are naturally lead to define a corresponding 
regularity notion intrinsically in Colombeau algebras as follows. 

Definition 13. Let 7 be an admissible scaling function and s be a real number. 
A Colombeau function U G QiW"^) is said to be globally of "f- Zygmund regularity 
s if for all a G there is M e Nq such that for aU (j) e ^a/(M") we can find 
positive constants C and 77 such that 

(15) |9"C/(0„x)|< xGR'",0<e<ry. 

I C 7e it |Q!| > s 

The set of all (globally) 7-Zygmund regular Colombeau functions of order s will 
be denoted by gt^^(W^). 

A detailed analysis of Ql ^ in arbitrary space dimensions and not necessarily 
positive regularity s will appear elsewhere. Here, as an illustration, we briefly 
study the case to = 1 and s > in some detail. Concerning applications to 
PDEs this would mean that we are allowing for media of typical fractal nature 
varying continuously in one space dimension. For example one may think of a 
coefficient function / in CJ(M) to appear in the following ways. 

Example 14. (i) Let / be constant outside some interval {—K, K) and equal 
to a typical trajectory of Brownian motion in [—K,K]] it is well-known 
that with probability 1 those trajectories are in C^R) whenever s < 
1/2. This is proved, e.g., in Sect. 4.4, elegantly by wavelet transform 
methods. 
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(ii) We refer to Sect. V.3, for notions and notation in this example. Then 
similarly to the above one can set / = in (—00, 0], / = 1 in [27r, 00) and 
in [0, 2tt] let / be Lebesgue's singular function associated with a Cantor- 
type set of order d £ N with (constant) dissection ratio < ^ < 1/2. Then 
/ belongs to C|(M) with s = log(d + l)/|log(^)|. (The classical triadic 
Cantor set corresponds to the case d = 2 and ^ = 1/3.) 



We have already seen that the Colombeau embedding does not change the mi- 
crolocal structure (i.e., the 7- wave front set) of the original distribution. We 
will show now that also the refined Zygmund regularity information is accurately 
preserved. If n S Np we denote by C^(R) the set of all n times continuously 
differentiable functions with the derivatives up to order n bounded. Note that 
C^"(]R) is a strict superset of C^'+^K)- 

Theorem 15. Let x G -4o(R) and s > 0. Define n G Nq such that n < s < n + 1 
then we have 



(16) ^lic^m n Gi.m = q{cm))- 

In other words, in case < s < 1 we can precisely identify those Colombeau 
functions that arise from the Zygmund class of order s within all embedded 
bounded continuous functions. 



Proof. We use the characterizations in [Q, Thms. 2.9.1 and 2.9.2 and the re- 
marks on p. 48 following those; choosing a smooth compactly supported wavelet 
g of order n we may therefore state the followingj^: / G (M) belongs to (M) 
if and only if there is C > such that 

(17) \Wgf{x, r)\ < Cr' for aU x. 



Now the proof is straightforward. First let / G C^(R). If |a| — 1 < n then 
=JfJ(5"/)l < II<9"/|Il~IIxIIli by Young's inequality. If \a\ > n we 
use lemma |lj and set 7^ = 7(eZ(0)) to obtain 



where Xa is a wavelet of order at least n. Hence (y_7|) gives an upper bound 
(77!"' * uniformly in x. Hence (15) follows. 

Finally, if we know that / G C^(K) and tj^(/) G SJ^^(R) then combination of 
( p^ and lemma |l^ gives if |a| > s 

l7l"IW^x»/(^,77')l = |a"4(/)(^e,x)| < C7i"l-^ 



uniformly in x. Hence another application of (17) proves the assertion. □ 



'^Note that we do not use the wavelet scahng convention adapted to L^-spaces here 
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